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Abstract — Unique word - orthogonal frequency division mul- 
tiplexing (UW-OFDM) is a novel OFDM signaling concept, 
where the guard interval is built of a deterministic sequence 
- the so-called unique word - instead of the conventional 
random cyclic prefix. In contrast to previous attempts with 
deterministic sequences in the guard interval the addressed UW- 
OFDM signaling approach introduces correlations between the 
subcarrier symbols, which can be exploited by the receiver 
in order to improve the bit error ratio performance. In this 
paper we develop several linear data estimators specifically 
designed for UW-OFDM, some based on classical and some 
based on Bayesian estimation theory. Furthermore, we derive 
complexity optimized versions of these estimators, and we study 
their individual complex multiplication count in detail. Finally, 
we evaluate the estimators' performance for the additive white 
Gaussian noise channel as well as for selected indoor multipath 
channel scenarios. 



I. Introduction 

In HI, El we introduced an OFDM signaling scheme, where 
the usual cyclic prefixes (CP) [3] are replaced by deterministic 
sequences, that we call unique words (UW). A related but - 
when regarded in detail - also very different scheme is known 
symbol padded (KSP)-OFDM E)-||6). Fig. QiHTc] compare the 
CP-, KSP-, and UW-based OFDM transmit data structures. 

In CP- as well as in UW-OFDM the linear convolution 
of the transmit signal with the channel impulse response is 
transformed into a cyclic convolution. Note that apart from 
the very first UW in the symbol stteam (see Fig. QJ) each 
UW plays a double role: The (i + l) th UW represents the tail 
of the i xh OFDM symbol, while it additionally represents the 
'cyclic prefix' for the (i + l) th OFDM symbol. However, there 
are some fundamental differences between the CP-based and 
the UW-based approach: 

> Different to the CP, the UW is part of the discrete Fourier 
transform (DFT)-interval as indicated in Fig. Q] Due to 
that reason the bandwidth efficiencies of UW-OFDM and 
conventional CP-OFDM are almost identical. 
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(a) Data structure using CPs 
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(c) Data structure using UWs 

Fig. 1: OFDM transmit data structures. 



> The CP is a random sequence, whereas the UW is 
deterministic. Thus, the UW can optimally be designed 
for particular needs like synchronization and/or channel 
estimation purposes at the receiver side. 

The broadly known KSP-OFDM uses a structure similar 
to UW-OFDM, since the known symbol (KS) sequence is 
deterministic as well. The most important difference between 
KSP- and UW-OFDM is the fact, that the UW is part of 
the DFT interval, whereas the KS is not. On the one hand 
this characteristic of the UW implies the cyclic convolution 
property addressed above, and on the other hand, but least 
that important, the insertion of the UW within the DFT in- 
terval requires to introduce some correlations in the frequency 
domain, which can advantageously be exploited by the receiver 
to improve the bit error ratio (BER) performance. Whilst 
in both schemes the deterministic sequences can be used 
for synchronization and channel estimation purposes, KSP- 
OFDM does not feature these correlations. We notice that 
KSP-OFDM coincides with zero padded (ZP)-OFDM Q if 
the KS sequence is set to zero. 

For single carrier/frequency domain equalization (SC/FDE) 
systems |8|-|22|, the benefits of UW based transmission have 
already sufficiently been studied 03]-|[l8), E0)-||2T]. The 
introduction of UWs in SC/FDE systems is straightforward, 
since the data symbols as well as the UW symbols are defined 
in time domain. In UW-OFDM the data symbols are defined 
in frequency domain, whereas the UW symbols are defined 
in time domain, which leads to some difficulties. In ll23l we 
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compared the similarities and differences of the UW approach 
for OFDM and SC/FDE. 

In our concept described in (T) we suggested to generate 
UW-OFDM symbols by appropriately loading so-called redun- 
dant subcarriers. The minimization of the energy contribution 
of the redundant subcarriers turned out to be a challenge. 
We solved the problem by generating a zero UW in a first 
step, and by adding the desired UW in a separate second step. 
We showed that this approach generates OFDM symbols with 
much less redundant energy [2] than a single step or direct UW 
generation approach as e.g. described in fl6l . Additionally, we 
optimized the positions of the redundant subcarriers to further 
reduce their energy contribution. We notice, that the concept in 
lTT6l generates completely different OFDM symbols compared 
to our approach in (TJ, and it has to deal with extremely 
high symbol energies and with the fact, that the performance 
depends on the particular shape of the UW. This is clearly 
in contrast to our approach, where the BER performance is 
independent of the particular shape of the UW due to the two- 
step generation approach. The BER behavior only depends on 
the freely selectable UW energy. 

The generation of the zero UW introduces a systematic 
complex valued block code structure within the sequence 
of subcarriers. From this point of view, the gain due to 
the exploitation of correlations in frequency domain can be 
regarded as a coding gain. Although it seems obvious at first 
glance to use an algebraic decoding approach, this decoding 
method fails due to the ill-conditioned nature of the linear 
system of equations to be solved l24l . ||251 . Instead, we 
showed that together with a linear minimum mean square 
error (LMMSE) data estimator ('decoder'), the concept shows 
a remarkable BER performance, particularly in frequency 
selective channels, where it clearly outperforms CP-OFDM 
lHJ. The performance can even be increased by allowing 
some systematic noise in the guard interval |26l . Several 
other attempts of applying UWs in OFDM systems can be 
found in the literature, e.g. in ll27l - ll28l . However, in all these 
approaches the guard interval and thus the UW is not part of 
the DFT interval. Therefore, in contrast to our UW-OFDM 
concept described below, no coding is introduced by these 
approaches. 

The aim of this paper is to give a comprehensive view 
on optimum and suboptimum linear data estimation principles 
particularly designed and optimized for UW-OFDM. We clas- 
sify the estimators into classical unbiased estimators and linear 
Bayesian estimators, respectively. We particularly investigate 
the theory of the estimators, and we give a comparison in 
terms of BER performance and in terms of a detailed study of 
the computational complexities. Furthermore, we emphasize 
the differences of the derived estimators to their counterparts 
in competing block oriented approaches like CP-OFDM and 
SC/FDE. The paper is organized as follows: In Sec. [TT] we 
briefly review the procedure of the unique word generation 
and the overall system model which has already been adressed 
in (TJ, J3, ll23l . and ll26ll . Next we derive data estimators 
for UW-OFDM using classical estimation theory approaches 
in Sec. [HI] leading to zero forcing (ZF) equalizer concepts. 



In particular we investigate the best linear unbiased estimator 
(BLUE) which represents a well known concept with a huge 
number of applications in engineering. In this work it is 
applied to UW-OFDM for the first time, and it turns out that its 
construction significantly differs from its counterparts in CP- 
OFDM and SC/FDE. Different to the BLUE for CP-OFDM 
and CP-SC/FDE, the determination of the BLUE requires the 
inversion of a full instead of a diagonal matrix. However, we 
derive a functionally equivalent but highly complexity reduced 
version of the BLUE within Sec. [Till In contrast to CP-OFDM 
where the BLUE represents the unambiguous zero forcing 
(ZF) solution, an infinite number of ZF solutions exists for 
UW-OFDM. We introduce two suboptimum low complexity 
ZF data estimators, that is the obvious channel inversion (CI) 
estimator and a quite intuitive estimator that we call time 
domain windowing (TDW) equalizer. Then, in Sec. IIVI linear 
Bayesian MMSE estimators are regarded. The basic version 
can already be found in (T). Similar as for the BLUE we 
derive a complexity reduced batch solution, and in addition we 
introduce a highly complexity optimized sequential version of 
the LMMSE estimator. It turns out that the latter features the 
lowest complexity of all regarded LMMSE estimator versions. 
In Sec. [V] we determine and compare the computational 
complexity of all presented data estimators. Finally, in Sec. I VII 
we highlight the BER performance of the introduced methods 
in the AWGN channel and in frequency selective indoor 
multipath environments. We conclude our work in Sec. IVHI 

Notation: Lower-case bold face variables (a, b,. . . ) indicate 
vectors, and upper-case bold face variables (A, B,. . . ) indicate 
matrices. To distinguish between time and frequency domain 
variables, we use a tilde to express frequency domain vectors 
and matrices (a, A,. . . ), respectively. We further use K to 
denote the set of real numbers, C to denote the set of complex 
numbers, I to denote the identity matrix, ( ) T to denote 
transposition, (-) H to denote conjugate transposition, E[-) to 
denote expectation, and tr(-) to denote the trace operator. For 
all signals and systems the usual equivalent complex baseband 
representation is applied. 

II. Review of UW-OFDM: Unique Word Generation 
and System Model 

We briefly review our approach of introducing unique words 
in OFDM time domain symbols, for further details see (TJ, J2). 
A block diagram of the transceiver chain is given in Fig. [2] 
Let x„ G C Ar " xl be a predefined sequence which we call 
unique word. This unique word shall form the tail of each 
OFDM time domain symbol vector of total length N. Hence, 
an OFDM time domain symbol vector, as the result of a length- 
A^-IDFT (inverse DFT), consists of two parts and is of the form 
[xj x^] , whereas only x^ € Cf- N ~ Nvr ' xl is random and 
affected by the data. In the concept suggested in (TJ, (2) we 
generate an OFDM symbol x = [xj T ] with a zero UW 
in a first step, and we determine the final transmit symbol x' = 
x + [0 T x^T] by adding the desired UW in time domain 
in a second step. As in conventional OFDM, the quadrature 
amplitude modulation (QAM) data symbols (denoted by the 
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Fig. 2: Block diagram of the UW-OFDM transceiver system. 



vector d E C NdXl ) and the N z zero subcarriers (at the band 
edges and at DC) are specified in frequency domain as part of 
the vector x, but here in addition the zero-word is specified 
in time domain as part of the vector x = F^- x. Here, Fjv 
denotes the length- 7V-DFT matrix with elements [Fjv]fcz = 
e -Jif k« f or k,l = 0, 1, N — 1. The system of equations x = 
F^ x with the introduced features can be fulfilled by spending 
a set of redundant subcarriers. We let the redundant subcarrier 
symbols form the vector r 6 C WrXl with N r = N u , we further 
introduce a permutation matrix P G C^ Nd+Nr ^ x ^ Nd+Nr \ and 
form an OFDM symbol (containing N z = N — Nd — N r zero 
subcarriers) in frequency domain by 



x = BP 



(1) 



B G (£Nx(N d +N r ) mser ( S me zero subcarrier symbols, and 
consists of zero-rows at the positions of the zero subcarriers, 
and of appropriate unit row vectors at the positions of data 
and redundant subcarriers. We will detail the reason for the 
introduction of the permutation matrix P and its specific 
construction shortly below. The time - frequency relation 



x can now be written as F^BP [d T r T ] 



J 



With 



M = F^BP = 



Mn M12* 
M 2 i M 22 



(2) 



where Mjy are appropriate sized sub-matrices, it follows that 



M 2 id 
matrix 



M 22 r = 0, and hence r = -M 22 M 2 id. With the 

T = -M 2 - 2 1 M 2 i G C Nr * Nd , (3) 



the vector of redundant subcarrier symbols can thus be deter- 
mined by the linear mapping 



Td, 



(4) 



cf. Fig. |2] The construction of T and thus also the energy 
of the redundant subcarrier symbols highly depend on the 
choice of P. The mean symbol energy E x > — E[x.' H x'] can 
be calculated to 

= ^(N d a 2 d + a 2 MTT H )) + xfx„, (5) 



Ex 



Er 



cf. 12. -S- and j£ describe the contributions of the data and 
the redundant subcarrier symbols to the total mean symbol 



energy before the addition of the UW, respectively, and E Xii 
describes the contribution of the UW. In U we suggested to 
choose P by a minimization of the symbol energy E x i which 
leads to the optimization problem 



argmin{tr(TT ff )} 



(6) 



where T is derived from (f3j) and ©, respectively. In Sec. 
[VTl we give an example of the optimum redundant subcarrier 
distribution for a specific parameter setup. 

With (O the vector c s = [d T f T ] of data and redundant 
subcarrier symbols can be written in the form 



c, = 



d = G,d. 



In (|7]l the matrix 



c (N d + N r )xN d 



(7) 



(8) 



can be interpreted as the code generator matrix for a systematic 
complex valued block code, that generates the code words c s . 
Note, that the subscript V in the code words c s and in the 
code generator matrix G s denotes sorted. In contrast to fl] 
and [2 | we do not incorporate the permutation matrix P into 
the code generator matrix in this work. As we will see later on 
this is essential for deriving low complexity receiver concepts. 
With (O and with the frequency domain version of the UW 
x M = Fjv [0 T xj] the transmit symbol can now also be 
written as 

x' = F Jv 1 (BPG s d + i u ). (9) 

After transmission over a dispersive (e.g. multipath) channel 
a received OFDM time domain symbol can be modeled as 



H c x' + n (10) 
= H c F w 1 (BPG s d + x u )+n, (11) 

where n G C Nxl represents a zero-mean Gaussian (time 
domain) noise vector with the covariance matrix cr^I, and 
H c G C NxN denotes a cyclic convolution matrix originating 
from the zero-padded vector of channel impulse response 
coefficients h c G C Nxl . After applying a DFT to obtain 
y r = FAry r , we exclude the zero subcarriers from further 
operation, which leads to the down-sized vector yd — B T y,. 
with y<j G (£(Nd+N r )xl. 



yd 



B T F w H c F- 1 (BPG s d + x u )+B T F 7V n. (12) 
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The matrix H c = FatHcF^ 1 is diagonal and contains the 
sampled channel frequency response on its main diagonal. 
H d = B t Ft V H c F^ 1 B with H d G c( N «+ N r)x( N «+ N r) 
is a down-sized version of the latter excluding the entries 
corresponding to the zero subcarriers. The received symbol 
can now be written in the form of the affine model 



Yd 



H d PG s d + H d B T x u 



B T FArn. 



(13) 



Note that (assuming that the channel matrix or at least 
an estimate of it is available) HdB T x„ represents the known 
portion contained in the received vector yd originating from 
the UW. As a first preparatory step we therefore subtract the 
UW influence to obtain the corrected symbol in the form of 
the linear model 



y c = Yd- H d B T x„ 
= H d PG s d + w, 



(14) 
(15) 



with the noise vector w = B T Fjvn. For the low complexity 
versions of the BLUE and the LMMSE estimator to be derived 
in the subsequent sections it additionally turns out to be quite 
advantageous to re-sort the receive vector by applying P T 
to separate the data subcarrier symbols and the redundant 
subcarrier symbols. The re-sorted receive vector y follows to 

y = P T y c (16) 
= P T H d PG s d + P T B T Fjvn. (17) 

With the re-sorted (and still diagonal) channel matrix H s = 
P T H C ;P and the noise vector v = P T B T Fjvn we finally 
arrive at the linear model 



y = H s G s d + v. 



(18) 



III. Classical Data Estimators - Zero Forcing 
Solutions 

In this section we consider classical unbiased data estima- 
tors of the form 

d = Ey, (19) 

where E e £N d x{N d +N, describes the equalizer. Note that 
in classical estimation the data vector is assumed to be 
deterministic but unknown. In order for the estimator to be 
unbiased we require 

E[d] = E[Ey] = EH s G s d = d. (20) 

Consequently, the unbiased constraint takes on the form 

EU e G s = I, (21) 

which is equivalent to the ZF criterion for linear equalizers. 
The solution to (fJTJi is ambiguous. To show this we consider 
a singular value decomposition of H S G S G N d+N r )xN d as 



H,G S = U 



V 



H 



(22) 



with unitary matrices U G c (Ar ' i+Ar '') )< ( Ard+Ar '') and V G 
C NdXNd , and with the diagonal matrix £ G R NdXNd having 



as its main diagonal the singular values of H S G S . With 
the unbiased constraint (or ZF criterion) (fJTJ becomes 



EU 



tH 



I 



(23) 



It is easy to see that (1231 and therefore also (fJTJ is fulfilled 
by every equalizer of the form 



E = V [5T 1 A] \J H 



(24) 



with arbitrary A G C NdXNr . We notice that the fact that 
the ZF solution is ambiguous distinguishes UW-OFDM from 
competing block oriented single input single output (SISO) 
approaches like e.g. CP-OFDM and CP-SC/FDE. For CP- 
OFDM the channel inversion receiver E = H^ 1 represents 
the unambiguous ZF solution which also corresponds to the 
optimum data estimator, cf. (3). For CP-SC/FDE the ZF 
solution is also unambiguous, it is given by the inverse of 
the diagonal symbol spaced channel matrix which contains 
the influence of the transmit pulse shaping filter, the dispersive 
(e.g. multipath) channel, and the receiver filter (e.g. a matched 
filter), cf. Ifl9l. 

Since the solution to the unbiased constraint is not unam- 
biguous it makes sense to look for the optimum solution which 
is commonly known as the best linear unbiased estimator. 

A. Best Linear Unbiased Estimator (BLUE) 

By applying the Gauss-Markov theorem |29l to (II St . and 



with the noise covariance matrix = E [vv H ] = Ncr^I, 
the BLUE and consequently the optimum ZF equalizer follows 
to 

Eblue - (Gf Hf H.G.J-^f Hf . (25) 

Eblue as given in d25T l represents the pseudoinverse of H S G S . 
Since the noise in (fT8l is assumed to be Gaussian, (|25T > 
is also the minimum variance unbiased (MVU) estimator. 

The covariance matrix of d = EeLUEy, or equivalently the 

covariance matrix of the error e = d d is given by 

C e - g - Nal{G* Hf HsGs)- 1 . (26) 

With the singular value decomposition as in (l22"i i. and after 
some rearrangements using standard matrix algebra, d25l > can 
immediately be re-written as 



E 



BLUE 



= V[£- 1 0]XJ H . 



(27) 



By comparing this result with (124b it can be concluded that 
Eblue corresponds to the solution in d24T > for the particular 
case A = 0. Eblue is in general a full matrix, which is in 
contrast to CP-OFDM and CP-SC/FDE, where the BLUE is 
given by a diagonal matrix. 

B. Complexity Optimized Version of the BLUE 

One drawback of the BLUE represented as in d25l l is the 
fact, that an Nd X Nd matrix has to be inverted to determine 
the equalizer. In this section we derive a significantly com- 
plexity reduced version of the BLUE by exploiting the simple 
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structures of G s 
decompose H s a 



and H^, respectively. For this purpose we 



H, 









H s .2 



(28) 



with the diagonal matrices H S) i e £ N d><N d an( j jj s2 
C NrXNr . With ® it follows that 



H„G S 



H M 




I 




H s .i 


H s2 




T 




H s2 T 



(29) 



and the expression (G^H^H S G S ) 1 appearing in d25l l and 
can be written as 



(30) 



(GfHfH.G,)" 1 = (H^H^+T^H^H^T)- 
We introduce the real diagonal matrices 



H s 2 H s 2 G 



^N r xN r 



(31) 
(32) 



and apply the matrix inversion lemma, cf. 11291 , to the right 
hand side of (f30b to obtain 

(Gf Hf H.G,)- 1 = 

D^ 1 - B^T" {TB^T 11 + DJ 1 )- 1 TD^ 1 . (33) 

The inversions of the real diagonal matrices Di and D 2 are 
trivial, and the additional matrix (TD^T ff + D^ 1 ) to be 
inverted is Hermitian and only has the dimension N r xN r . Fur- 
thermore, the expression TD^ 1 (and its Hermitian transpose) 
occurs repeatedly in d33l which allows for further complexity 
reduction. 

In section [V] we will study the complexity of the different 
representations of the BLUE. We note that the derivation of 
the complexity reduced version of the BLUE has mainly been 
made possible by the re-sorting (multiplication with P T ) of 
the data and redundant subcarrier symbols in (fTTT i. 



C. Sub-Optimum ZF Receiver Structures 

Any unbiased linear data estimator, or equivalently, any 
linear zero forcing equalizer has to fulfill (fJTJ. As already 
shown above the ZF solution is ambiguous for the UW-OFDM 
transmission model described in ( fT8l . Another quite intuitive 
and straightforward ZF solution is given by 



Eci = [I o] h; 1 . 



(34) 



This equalizer inverts the channel H s first, and the data 
symbols are extracted subsequently. Clearly this procedure 
fulfills (f2Tb . In the following we will refer to this equalizer as 
the channel inversion (CI) receiver. Using the decomposition 
of H s as in ((28), d34t can be simplified to 



Eci = (H Sj i) _1 [I 0] 



(35) 



The channel inversion receiver represents a low complex 
solution since H S) i has a diagonal structure, but it does not 
take advantage of the correlations introduced by G s at the 
transmitter side. The covariance matrix of d = Eciy> or 



equivalently the covariance matrix of the error e 
can easily shown to be 

Cee = ^(HfiH.,!)- 1 . 



(36) 



Next we address another quite intuitive equalizer that ex- 
ploits the a-priori knowledge, that the guard interval samples 
of an UW-OFDM symbol must be zero after the channel 
inversion in the noiseless case. In the presence of noise we 
therefore simply force the guard interval samples to zero which 
is achieved by an equalizer of the form 



E. 



TDW 



[I 0]P T B T F 



where 



W 



N 



I 




WF^BPH; 1 , 



(37) 



(38) 



The time domain windowing (TDW) equalizer starts with an 
inversion of the channel, next the permutation is applied and 
the zero subcarrier symbols are added again in order to be able 
to transform back to time domain with a length- 7V-IDFT. Here 
a windowing (described by W) takes place, where the guard 
interval samples are forced to zero. Next a transformation 
back to frequency domain is performed, the zero subcarriers 
are excluded again, a re-sorting is done, and finally the data 
symbols are extracted. It can easily be shown, that Etdw 
also fulfills (|2TT >. Note that the TDW equalizer also represents 
a quite low complex solution since none of the individual 
operations requires a full matrix multiplication, in fact most of 
the steps apart from DFT and IDFT are trivial. The covariance 

matrix of d = EtdwY. or equivalently the covariance matrix 
of the error e = d — d is given by 



Cee = Na^ Etdw E?dw 



(39) 



IV. Linear Bayesian Data Estimators - LMMSE 
Solutions 

We now turn to the widely used linear minimum mean 
square error data estimator which is derived with the help 
of the Bayesian approach. In the Bayesian approach the data 
vector is assumed to be the realization of a random vector 
instead of a deterministic and unknown vector as in the 
classical estimation theory applied above. In the following 
we derive the LMMSE batch solution, next we formulate a 
complexity optimized version of the LMMSE batch solution, 
and finally we derive a highly complexity optimized version 
of the sequential LMMSE estimator. 



A. LMMSE Batch Solution 

By applying the Bayesian Gauss-Markov theorem [29 1 to 
dl8t . where we now assume d to be the realization of a 
random vector, and by using Cjj = (TjI and C^c = Na^l 
the LMMSE equalizer follows to 



E 



LMMSE 



wh: 



(40) 
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where W represents a Wiener smoothing matrixQ given by 



W = G? I G S G 



iH 



(41) 



(|40T > allows the following interpretation of the LMMSE esti- 
mator's mode of operation: The LMMSE equalizer acts as a 
composition of a simple channel inversion stage (multiplica- 
tion with H.T 1 as in d34l i) and a Wiener smoothing operation 
(multiplication with W). The Wiener smoothing operation 
exploits the correlations between subcarrier symbols which 
have been introduced by (0]i at the transmitter, and acts as a 
noise reduction operation on the subcarriers. For the equalizer 
in (00]!, an (N d +N r ) x (N d +N r ) matrix has to be inverted. By 
applying the matrix inversion lemma, it can easily be shown 
that the equalizer can equivalently be determined by 

Elmmse = (G^H^H S G S H J^I) 1 G^H^ . (42) 

d42l shows strong similarities to the BLUE in d25b . For <r^ = 
the expressions for the LMMSE equalizer and the BLUE 
coincide. Note that by using (l42l instead of d40b for the 
LMMSE equalizer determination the matrix to be inverted only 
has the dimension Nd x N d . The error e = d d has zero 
mean, and its covariance matrix is given by 



JV^( G fHfH s G 



-I) 



(43) 



B. Complexity Optimized LMMSE Batch Equalizer 

For the LMMSE equalizer a complexity reduced version 
can be derived similar as for the BLUE in Sec. IIII-BI By 
introducing the real diagonal matrices 



Di = HfjHs.i H =2-1, 

D2 = H^ 2 H S ,2, 



(44) 
(45) 



the expression (G^H^H S G S + ^p-T) 1 appearing in d42l 
and in d43l can be written as 

(G?HfH.G a + ^I)- 1 = 

D" 1 - D^T^ (TD^T^ + 1 ) -1 TD^ 1 . (46) 
The derivation widely coincides with the one in Sec. IIII-BI 



C. Complexity Optimized Sequential LMMSE Receiver 

In this section we derive a highly complexity optimized 
sequential LMMSE receiver. We address the equalizer deter- 
mination procedure as well as the data estimation process. 
The sequential LMMSE estimator completely avoids matrix 
inversions. Again the preparatory steps described in section 
HIl - especially the re-sorting in (fTTT i - are extremely beneficial 
for the derivation of the complexity optimized solution. In this 
section we use the system model 



y = H s c s 



v. 



(47) 



'Even though we use the same notation the Wiener smoothing matrix has 
nothing to do with the matrix in i38t . 



and estimate c s which includes both the data and the redundant 
subcarrier symbols. It turns out that using the system model 
in d47b instead of the one in ( fTST l drastically simplifies the 
sequential LMMSE procedure since H s is diagonal, in contrast 
to H S G S . We let c s [n] be the LMMSE estimate based on the 
first n + 1 elements {y[0], . . . , y[n]} of the vector y, and 
M[n] be the corresponding minimum MSE matrix 



M[n] = E[(c s -5>])(c s -Z s [n}) H }. 



(48) 



Furthermore, h s [n] denotes the column vector that corresponds 
to the Hermitian transpose of the n th row of H s [n]. The 
sequential LMMSE estimator for the Bayesian linear model 
as in (07]) becomes (cf. Il29l ): 



Initialization: 

c,[-l] 
M[-l] 



E[c a ] = 

Gc„c s 

= E \{c s — c s [— l])(c s - E s [-l]) ff 
I T H 



Od^*s^s = a d 



(49) 



(50) 



For n = 0, 1, 



Gain Vector Update: 



(N d + N r -l) do 

M[n-l]h s [n] 



o? + hf [n]M[n- l]h s [n] 



(51) 



Minimum MSE Matrix Update: 



M[n] = (I - k[n]hf [n])M[n - 1} 



(52) 



Estimate Update: 



c. n = cJn - 1 



k[n](y[n]-h?[n]c.[n-l]) (53) 



( T5Tb and d52b can be regarded as the equalizer determination 
procedure that can completely be performed immediately after 
channel estimation. Note that only the final MSE matrix 
M[N d +N r -l] but all gain vectors (k[0], k[l], . . . , k[N d + 
N r — 1]) are required to be stored until the next channel 
estimation update. d53l describes the sequential data estimation 
procedure for one UW-OFDM symbol, that has to be applied 
to every received OFDM symbol. After (N d + N r ) iterations 
the vector ^ 

5=[I 0]6 S , (54) 

that contains the first N d entries of c s , exactly corresponds to 
the data estimate obtained when applying the batch LMMSE 
equalizers ( |40T > or d42b . Further, the upper left N d x Nd sub- 
matrix of M[N d + N r - 1] 



Cgg - [I 0] M[N d + N r - 1] 



(55) 



exactly corresponds to the error covariance matrix in d43l ). 
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In the following we significantly simplify the sequential 
LMMSE procedure. We first exploit the fact that the system 
matrix H s is diagonal. Let [H s ] n „ be the n th main diagonal 
element of H s , [M] nn be the n th main diagonal element 
of M, rh„ be the n th column of M, and c s _ n be the n th 
element of c s . Then since H s is diagonal the iteration steps 
can be simplified to: 



Gain Vector Update: 



k[n] = 



[H s ]* nn m n [n - 1] 



<r?, + \m a 



! [M] nn [n-l] 



(56) 



Minimum MSE Matrix Update: 

M[n] - M[n - 1] - [H s ] nn k[n]xn%[n - 1] (57) 

Estimate Update: 

c s [n] = c s [n - 1] + k[n](y[n] - [H s ] n „ ■ c s , n [n - 1]) (58) 

We notice that in the MSE matrix update equation a full matrix 
multiplication simplifies to a (column x row) multiplication, in 
the gain vector update equation a matrix-vector multiplication 
simplifies to a vector-scalar multiplication, and in the estimate 
update equation a vector inner product simplifies to a scalar 
product. 

The equations for the first Nd iteration steps can further 
significantly be simplified by exploiting the fact that the data 
symbols are mutually uncorrelated, i.e. the upper left Nd x 
Nd submatrix of M[— 1] is diagonal (and real), cf. d50b . We 
partition the gain vector and the MSE matrix as 



k = 



;M = 



M dr 



(59) 



where the indices V and V indicate, that the corresponding 
vector and matrix entries correspond to data and redundant 
subcarrier symbols, respectively. Furthermore, we split all 
update equations in separate equations for the data and the 
redundant subcarrier symbols. The following consequences of 
Md[— 1] = ctjI can be exploited: 

• For all n = 0, 1, . . . , (Nd — 1) the gain vector kd[n] 
is non-zero only at its n th entry, and can therefore be 
replaced by the scalar gain factor kd[n]. Consequently, 

only one data symbol d n will be updated at the n th 
iteration step, so the estimate update equation for the 
data entries simplifies to a scalar equation. Since the data 

estimates are initialized with zeros (d n [— 1] = 0), the data 
estimate update equation becomes particularly simple. 
« For all n = 0, 1, . . . , (Nd — 1) the matrix M^[n] is 
diagonal and real, and at the n th iteration step (again 
for n — 0,1, ... , (Nd — 1)) it only needs to be updated at 



its n main diagonal element [MJ 



Consequently, 



the MSE matrix update equation for M^[n] also reduces 
to a scalar equation. 



• Due to similar arguments the matrix update for M^n] 
simplifies to an update of its n th column vector rhdr,n [n] 
for 7i = 0, 1, ... , (Nd — 1). 
The iteration equations for n = 0,1, (Nd — 1) finally 
simplify as follows: 

Gain Vector Update: 



kd\n] 



k r [ii] = 



[H s ];„rh dr ^[n - 1] 

o2 + o3i[fi a ] nn p 

Minimum MSE Matrix Update: 

[M d ] nn [ri] = 



(60) 
(61) 



fH, 



(62) 



(63) 



mdr,n[n] = m drtn [n - 1] - o^[H s ]„„k r [n] 

M r [n] = M r [n - 1} - [U s } nn k r [n]mgjn - 1] (64) 



Estimate Update: 



d n = k d [n]y[n] 
r[n] = r[n — 1] + k r [n]y[n] 



(65) 
(66) 



For the first Nd iteration steps the highly complexity re- 
duced equations d60b to d64l ) can be used for the equalizer 
determination, only for the last N r iteration steps the more 
complex (but also quite simplified) equations d56l > and ( T57l ) 
have to be evaluated. Similarly the first Nd iteration steps of 
the data estimation procedure for an UW-OFDM symbol can 
be performed using the highly complexity reduced equations 
( l65b and d66l l. while for the last N r steps ( 158b has to be used. 
The complexity analysis will be given in section [V] 

The derived procedure also allows for a quite intuitive 
interpretation of the mode of operation of the sequential 
LMMSE estimator: For n = 0,1,..., (N d - 1) only one 
data symbol is updated in each iteration step. Consequently, 
during the first Nd iterations we only count one single complex 
multiplication per data subcarrier symbol as in classical CP- 
OFDM. Merely the redundant symbols (which require a vector 
update) are truly updated from step to step. Only for the 
last N r iteration steps also the data subcarrier symbols are 
updated from iteration to iteration by utilizing the correlation 
information contained in the redundant subcarrier symbols. 

Note that these simplifications would not have been possible 
without the re-sorting step in ( fPTI ). Without the re-sorting 
step the gain vector would already be filled completely in a 
very early iteration step, namely immediately after the first 
redundant subcarrier symbol appears within y. Furthermore, 
if we had used the system model (fT8l instead of d47l i. then 
we would have to perform full (Nd x Nd) • (Nd x Nd) matrix 
multiplication operations for the last 7Y r iteration steps of the 
MSE matrix update. 



s 



V. Complexity analysis 

In this section we will analyze the computational com- 
plexity of the derived equalizers on the one hand, and of 
the corresponding data estimation procedures on the other 
hand. These investigations will clearly show the benefits of the 
complexity reduced versions. In practice, the equalizers need 
to be determined each time the channel estimate is updated. 

A. Prerequisites 

We are aware that it is difficult or even impossible to declare 
an equitable measure of complexity, since the complexity of 
an implementation strongly depends on the choice of the 
hardware and software architecture and of many implemen- 
tation details. Some operations can even be implemented in 
many different ways, which might have advantages on certain 
architectures as well. To simplify things we basically count 
the number of complex multiplication equivalents (CME) for 
each individual equalizer and for the corresponding data esti- 
mation procedure. We completely ignore additions. Complex 
division are counted as 1 CME. Since the number of required 
divisions is negligible, this simplification does not effect the 
final complexity considerably. Real multiplications and real 
divisions are counted as \ CME. 

For many of the derived equalizer implementations we have 
to deal with matrix products of the form A _1 B with a positive 
definite Hermitian matrix A e C mxm and with B e C mx " 1 . 
We notice that calculating X = A _1 B is equivalent to solving 
the systems of simultaneous linear equations 

AX = B. (67) 

For our complexity calculations we assume that d67l ) is solved 
with the help of a Cholesky decomposition of A given by 
A = LL H , where L is a lower triangular matrix having 
positive values on its main diagonal. AX = B can be rewritten 
as L(L H X) = B. To obtain X one can solve LY = B for 
Y with the help of a forward substitution, and subsequently 
solve L ff X = Y for X with the help of a backward sub- 
stitution. The Cholesky decomposition requires |to 3 complex 
multiplications/divisions and m square roots [30|, [31 1. We 
neglect the square roots and end up with g m 3 CME. A single 
forward or backward substitution requires ^m 2 + \m CME. 
To solve ( |67T ) with the help of a Cholesky decomposition we 
can finally assume a total count of |m 3 + m 2 rib + mni, CME. 

Whenever possible, we take any simplifications into ac- 
count, that a special matrix structure (e.g. a diagonal, a real 
or a Hermitian matrix) could offer. Exemplarily, if the result 
of a matrix product is Hermitian, e.g. as in X = A H A, then 
only the main diagonal and the lower triangular part needs to 
be computed. 

B. Complexity of the Investigated Equalizers and Data Esti- 
mation Procedures 

Before performing the data estimation with the help of 
one of the investigated equalizers an OFDM symbol has to 
be transformed to frequency domain with a length- iV-FFT 
(fast Fourier transform) which requires ^N\og 2 {N) CME. 



Furthermore, as one of the preparatory steps the influence 
of the UW has to be subtracted as described in ( fl4] i. Since 
we do not count additions/subtractions in our complexity 
considerations this step does not increase the CME count for 
the data estimation procedure. 

In the following we consider the complexity of the equaliz- 
ers investigated above. We start our complexity investigations 
with the most simple equalizer Eci as given in (|35T >. To 
determine Eci only Nd CME (namely complex divisions to 
invert H s ,i) are required. The data estimation procedure for 
an OFDM symbol in frequency domain requires Nd CME 
(namely complex multiplications). 

To estimate the data part of an OFDM symbol with the help 
of Etdw one could first determine its matrix representation 
as in d37T i, and then estimate the data vector by performing 

the full matrix-vector product d = EtdwY which requires 
N d (N d + N r ) operations. However, most of the individual 
operations required to perform the data estimation are trivial. 
The procedure starts with the multiplication H^y (Nd + N r 
CME), next the permutation is applied and the zero subcarrier 
symbols are added (zero CME) in order to be able to trans- 
form back to time domain with a length- A^-IFFT log 2 (iV) 
CME). In time domain a windowing takes place, where the 
guard interval samples are forced to zero (zero CME). Next 
a transformation back to frequency domain is performed 
(f log 2 (iV) CME), the zero subcarriers are excluded again, 
a re-sorting is done, and finally the data symbols are extracted 
(zero CME). So in total the data estimation procedure per 
OFDM symbol requires 7Y log^N) +N d + N r CME. For the 
equalizer determination only H s needs to be inverted which 
requires N d + iV r CME (namely complex divisions). 

Next we investigate the complexity of the different BLUE 
and LMMSE estimator batch representations. For all imple- 
mentations the data vector estimation for one OFDM symbol 
requires a full matrix vector product d = Ey with Nd(Nd + 
N r ) CME. The complexity of the equalizer determination 
differs significantly for the different implementations. We start 
with the representation of the BLUE as in ( f25l l and with the 
LMMSE estimator as in (l42l i. These two expressions merely 
differ in the regularization term which only adds a single 
arithmetic operation. We neglect this single operation and treat 
(|25| | and ( |42T > as equally complex. Using^ d29i > it is easy to 
see that the matrix multiplication Xi = H S G S only requires 
N d N r CME. For the product X 2 _= X^X X we can use the 
findings from namely X 2 = H^H S i +T H U^ 2 H S ^T. 
By additionally exploiting the fact that X2 is Hermitian, one 
can easily find that the matrix product X 2 = XfXi requires 
\NjN r + N d N r + N d + N r CME. Finally the operation 
(X 2 ) _1 Xf requires |iV 3 + N%N r + Nj + N d N r CME by 
using the Cholesky decomposition together with the forward 
and backward substitutions as mentioned above. The overall 
CME count for the BLUE in d25j and the LMMSE estimator 
in d42b therefore adds up to 

lN 3 d + ^N 2 d N r + 3N d N r + N 2 d + N d + N r CME. (68) 
With similar considerations one can show that the LMMSE 



9 



equalizer as expressed in (|40T > requires 

+ %N 2 d N r + 2N d N? + ±A/ r 3 

+ Nj + ^N d N r + |iV rf + |JV r CME. (69) 

For the complexity optimized batch representations of the 

BLUE and the LMMSE estimator one has to determine the 
expressions in d33l l and d46i l. respectively, followed by a matrix 
multiplication with G^H^. The simple inverses D^ 1 and 
D^ 1 require |7V d + |A^ r CME, to determine D^T^ another 
\N d N r , and for T\D^T H ] additional \N d N% + \N d N r 
CME are required. The operation (•) _1 TD^ 1 demands hN%+ 
N d Nt + N d N r CME, and the multiplication with D^T^ 
adds NjN r CME. The determination of Gf Hf and the final 
multiplication add N d N r + NjN r + Nj CME, cf. ([29), which 
totals to 

±A/ 3 + 2NjN r + lN d N? 

+ Nj + 3N d N r + § N d + f iV P CME. (70) 

Finally, we investigate the complexity of the sequential 
LMMSE estimator regarded in section HV-CI We start with the 
data vector estimation of one OFDM symbol. For the first N d 
iteration steps the estimate updates are performed using (l65l l 
and (l66l l which in total requires N d N r + N d CME (complex 
multiplications). For the last N r iteration steps (l58l has to 
be evaluated which in total requires N d N r + + N r CME 
(complex multiplications). Consequently the equalization of 
one OFDM symbol requires 

2N d N r + N? +N d + N r CME. (7 1 ) 

For the equalizer determination we have to count the oper- 
ations required for the gain factor updates and for the MSE 
matrix updates. For the first N d iteration steps the gain factor 
updates are performed using ( |60l and d6"T| i which in total 
requires N d N r + N d complex multiplications and 2N d real 
divisions (which we count as ^N d CME). For the last N r 
iteration steps ( f56l l has to be evaluated which in total requires 
N d N r + TV, 2 + N r complex multiplications, N r real multipli- 
cations (\N r CME) and 2N r real divisions (±N r CME). The 
MSE matrix updates for the first N d iterations are performed 
using ((62]i to (O which in total requires N d N? + N d N r 
complex multiplications and Nd real divisions (jN d CME). 
For the last N r iteration steps (l57T i has to be evaluated which 
in total requires NjN r + 2N d N? +N? + N d N r + TV 2 complex 
multiplications. We finally arrive at 

NjN r + 3N d N? + iV r 3 

+ 4N d N r + 2N? + lN d + lN r CME. (72) 

C. Numerical Example 

In the simulation section we will show results for a par- 
ticular parameter setup. The most important parameters can 
be found in Tab. HI] For the complexity considerations only 
N d ,N r and N are important. The particular choices in Sec. 
|Vl]are N d = 36, N r = 16, N = 64. Tab. U compares the 
complexity of the different equalizer representations and data 



estimation procedures, respectively, for that particular param- 
eter setup. Note that for the data vector estimation per OFDM 
symbol we count the contribution of the FFT (-ylog 2 (-/V) 
CME) which is required in all cases, and the additional effort 
contributed by the particular equalization procedure. 



TABLE I: Computational complexity of the introduced equal- 
izers and data estimators. 





CME for equalizer 


CME for data est. 


Equalization method 


determination 


per OFDM symbol 


Eci mi 


36 


228 


Etdw 02) 


52 


628 


Elmmse (!!!) 


127677 


2064 


Eblue,E LM mse J25), (42) 


88612 


2064 


Eblue,E L mmse (33). (46) 


59068 


2064 


Sequential LMMSE 


55387 


1652 



We observe, that the simple equalizers Eci and Etdw show 
a significantly lower complexity for the equalizer determina- 
tion as well as for the data estimation per OFDM symbol. 
Concerning the BLUE and the LMMSE estimator we can 
state that the complexity optimized batch solutions reduce 
the equalizer determination complexity by around 33% com- 
pared to the straightforward implementations in d2"5i l and d42l . 
respectively. The complexity optimized sequential LMMSE 
estimator which completely avoids matrix inversions further 
reduces the equalizer determination complexity by another 6%, 
and interestingly enough also the data estimation complexity 
can be reduced by 20%. 

VI. Simulation Results 

In this section we evaluate the introduced receiver concepts 
in terms of their BER performance. We notice that all derived 
variants of an estimator (e.g. of the LMMSE estimator) 
perform equivalently. Different performance of distinct ver- 
sions of an estimator would only be expected if fixed point 
implementations were regarded which is not the focus of our 
investigations. 

A. Simulation Setup 

We show simulation results with and without outer channel 
coding. For the case when an outer channel code is used, 
the block diagram in Fig. |2] is extended by an outer channel 
encoder and an interleaver at the transmitter side, and by a 
deinterleaver and decoder at the receiver side. We used the 
same outer convolutional encoder with the industry standard 
rate 1/2, constraint length 7 code with generator polynomials 
(133, 171) as defined in [32 1 . A soft decision Viterbi algorithm 
is applied for decoding. The main diagonal of the appropriate 
matrix Cgg is used to specify the varying noise variances along 
the data symbols after data estimation. We assumed perfect 
channel knowledge in the simulations to be presented below. 

In |fl] we compared our UW-OFDM approach with the 
CP-OFDM based IEEE 802.11a WLAN standard [32| and 
showed that UW-OFDM outperforms CP-OFDM in frequency 
selective indoor environments. In this work we use the same 
parameter setup as in (T| which has been adapted to the 
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TABLE II: Main PHY parameters of the investigated UW- 
OFDM system. 



Modulation scheme 


QPSK 


Coding rates (outer code) 


uncoded, 1/2 


FFT length 


64 


Occupied subcarriers 


52 


Data subcaniers 


36 


Redundant subcarriers 


16 


DFT period 


3.2 


Guard duration 


800 ns 


Total OFDM symbol duration 


3.2 (j,s 


Subcarrier spacing 


312.5 kHz 



802.11a standard wherever possible. The most important pa- 
rameters are specified in Table HI] The sampling frequency 
has been chosen to be f s = 20 MHz. As in ll32l the indices 
of the zero subcarriers within an OFDM symbol x are set to 
{0, 27, 28, ...,37}. The indices of the redundant subcarriers are 
chosen to be {2, 6, 10, 14, 17, 21, 24, 26, 38, 40, 43, 47, 
50, 54, 58, 62}. This set (which can also be expressed by an 
appropriate permutation matrix P) minimizes the cost function 
in (0, and therefore also the mean energy of the redundant 
subcarriers. Since we focus on data estimation procedures in 
this work rather than on synchronization or channel estimation 
approaches we chose the zero UW for the BER simulations 
below. Note that in conventional CP-OFDM like in the WLAN 
standard, the total length of an OFDM symbol is given by 
Tdft + To j. However, the guard interval is part of the DFT 
period in the UW-OFDM approach which leads to significantly 
shorter total symbol durations. Hence, the compared systems 
show almost identical bandwidth efficiencies. 

B. Simulation Results in the AWGN Channel 

Clearly, OFDM is designed for data transmission in fre- 
quency selective environments. Nevertheless, we start our 
comparison with simulation results in the AWGN channel, 
since these results provide first interesting insights. In Fig. [3] 
the BER performance of the different data estimators is 
compared under AWGN conditions. As in all following BER 
figures we present curves for the case no outer code is used 
(we label it 'uncoded' in the figures), and for an outer coding 
rate r = |. 

We start the discussion with the uncoded case: As ex- 
pected the CI estimator shows the worst performance, since it 
completely ignores the information present on the redundant 
subcarriers. Surprisingly, the very simple and intuitive TDW 
data estimator performs almost as well as the BLUE and 
the LMMSE in the AWGN environment. At a BER of 10~ 6 
these three estimators which all make use of the a-priori 
knowledge introduced by the zero UW outperform the CI 
estimator by around 1.5dB. The trend is similar for r = |. 
However, it is completely in contrast to single carrier systems 
(e.g. SC/FDE) that the LMMSE estimator and the BLUE 
show a different performance in an AWGN environment. This 
comes from the fact that in single carrier systems the received 
QAM symbols are uncorrelated in an AWGN environment, 




E b /N (dB) 

Fig. 3: Simulated BER performance of the investigated data 
estimators in the AWGN channel. 



whereas in UW-OFDM systems correlations are inherently 
present due to (0]). However, the performance gain of the 
LMMSE estimator is quite small, and the BLUE approaches 
the LMMSE estimator performance for high E^/Nq, as the 
term in ( f42l > converges to zero. 

C. Simulation Results in Frequency Selective Indoor Environ- 
ments 

For the simulation of indoor multipath channels we applied 
the model described in (33), which has also been used dur- 
ing the IEEE 802.11a standardization process. The channel 
impulse responses are modeled as tapped delay lines, each 
tap with uniformly distributed phase and Rayleigh distributed 
magnitude, and with power decaying exponentially. The model 
allows the choice of the channel delay spread. For a more 
detailed description we refer to ||33ll . For illustration purposes 
we use two different channel snapshots in this section, each 
channel featuring a delay spread of 100 ns, and a total duration 
not exceeding the guard interval. The frequency responses are 
shown in Fig. |4] Channel A does not show any deep fading 
holes, whereas channel B features two spectral notches within 
the system bandwidth, one at a data subcarrier position, the 
other one at a redundant subcarrier position. 

Let us first interpret the results for channel A, cf. Fig. 
We observe similar trends as in the AWGN case, but now the 
LMMSE estimator and the BLUE clearly outperform the TDW 
estimator. For uncoded transmission the TDW outperforms the 
CI estimator by 1.9dB (again at a BER of 10~ 6 ), the BLUE 
and the LMMSE estimator gain 2.6dB and 2.7dB, respectively. 
For r = 5 the corresponding gains shrink to l.OdB, 1.35dB 
and 1.65dB, respectively. 

Finally Fig. [6] shows the simulation results for channel B 
with its deep spectral notches. Very noticeable in the uncoded 
transmission is the bad performance of the CI and the TDW 
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Channel A 




Frequency (MHz) 



Fig. 4: Frequency responses of indoor multipath channel 
snapshots. 



o 




E b /N (dB) 



Fig. 5: Simulated BER performance of the investigated data 
estimators for channel A. 



estimators. Here the performance gain of the BLUE and the 
LMMSE estimator is significant. The performance of the CI 
estimator is dominated by the weak BER behavior of data 
subcarrier symbols corresponding to deep spectral notches in 
the channel frequency response, while the LMMSE estimator 
(and similarly the BLUE) considerably decrease the noise 
on that subcarriers. (They decrease the noise variance on all 
subcarriers, but the effect is significant on subcarriers corre- 
sponding to deep spectral notches, cf. [1]). The BLUE and the 
LMMSE estimator perform almost equivalently, which is again 
in contrast to SC/FDE systems, where the performance gain of 
the LMMSE estimator over the BLUE in channels with deep 
fading holes is usually much larger, particularly at low E b /No 
values, cf. lfT4l . In coded transmission the performance loss 
of the CI estimator compared to the best performing LMMSE 
estimator decreases to 1.7dB. The significant improvement 




E b /N (dB) 



Fig. 6: Simulated BER performance of the investigated data 
estimators for channel B. 



of the CI estimator in the coded case was expected as this 
corresponds to the usual coding gain as it is also observed 
in CP-OFDM. Somewhat unexpected, and in contrast to the 
uncoded results and those in an AWGN channel and in channel 
A, the TDW equalizer performs almost 0.7dB worse compared 
to the CI equalizer at a BER of 10~ 6 . To understand this effect 
we will now have a closer look on the way the TDW estimator 
works. In fact, although it is hardly noticeable in Fig. [6] in the 
uncoded case the TDW only outperforms the CI estimator in 
the high E^/Nq range, but performs worse in the low E^/Nq 
range (0-15dB). However, this is the interesting E^/Nq range 
for coded transmission. We will now have a look on the noise 
variances (after equalization) and later on the BERs on the 
individual data subcarriers. 

Fig. [7a] and [7b] show the normalized noise variances after 
equalization at a fixed E b /N (E b /N = 4 dB) for both data 
estimators. We observe that on the data subcarrier with index 
11 the noise variance is tremendously reduced by the TDW 
compared to the CI estimator. This data subcarrier corresponds 
to the deep spectral notch around 5 MHz in the channel's 
frequency response. However, we also notice that the noise 
variances on data subcarriers around data symbol No. 1 1 are 
a little bit higher for the TDW compared to the CI estimator. 
On average (when averaged over all data subcarriers) the TDW 
equalizer clearly reduces the noise power compared to the CI 
equalizer, but besides a significant noise reduction on highly 
attenuated subcarriers, the TDW equalizer 'distributes' some 
noise onto neighboring subcarriers. Fig. 1751 additionally shows 
the difference between the resulting BERs of the TDW and 
the CI estimators on a subcarrier basis. We observe, that the 
tremendous noise reduction by the TDW equalizer on the 
II th data subcarrier indeed leads to a lower subcarrier BER 
compared to the CI equalizer, but the improvement is minor. 
In return, the higher noise variances on the adjacent data 
subcarriers lead to increased corresponding subcarrier BERs 
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Fig. 7: Subchannel noise variances after CI and TDW data 
estimation, and BER difference per subchannel. 



for the TDW estimator. In total the increase of these subcarrier 
BERs lead to a worse overall BER performance of the TDW 
compared to the CI estimator for these E^/No values. The 
overall noise reduction by the TDW estimator is not translated 
to an overall BER gain for that particular channel for these 
Eb/N values. 

VII. Conclusion 

In this work we investigated several linear data estimators 
specifically designed for UW-OFDM. We introduced data 
estimators following the principles of classical estimation 
theory which lead to ZF equalizers. Two simple and intuitive 
ZF equalizers and the optimum ZF equalizer corresponding 
to the BLUE have been discussed. Following the Bayesian 
estimation principle the LMMSE estimator has been presented, 
and its batch and sequential versions have been regarded. We 
derived highly complexity reduced versions of the individual 
estimators and investigated their complexity in detail in terms 
of equivalent complex multilication counts. The CME count 
of the complexity optimized BLUE and LMMSE estimator 
versions could considerably be reduced compared to their 



straightforward counterparts, but still they show a significantly 
higher CME count compared to the simple ZF solutions. 
With the help of simulations we demonstrated the bit error 
behavior of the proposed estimators in the AWGN channel 
and in frequency selective indoor environments. Especially in 
frequency selective channels featuring deep fading holes the 
BLUE and in particular the LMMSE estimator significantly 
outperform the simple ZF estimators. 
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